arXiv:1507.07873v6 [physics.ao-ph] 6 Feb 2016 


Generated using the official AMS MbX template—two-column layout. FOR AUTHOR USE ONLY, NOT FOR SUBMISSION! 


Journal of Physical Oceanography 


Are rogue waves really unexpected? 

Francesco Fedele* 

School of Civil and Environmental Engineering, School of Electrical and Computer Engineering, Georgia Institute of Technology, 

Atlanta, GA, USA. 


ABSTRACT 


An unexpected wave is defined by |Gemmrich and Garrett| p008^ as a wave that is much taller than a set 
of neighboring waves. Their definition of ’’unexpected” refers to a wave that is not anticipated by a casual 
observer. Clearly, unexpected waves defined in this way are predictable in a statistical sense. They can occur 
relatively often with a small or moderate crest height, but large unexpected waves that are rogue are rare. Here, 
this concept is elaborated and statistically described based on a third-order nonlinear model. In particular, 
the conditional return period of an unexpected wave whose crest exceeds a given threshold is developed. 
This definition leads to greater return periods or on average less frequent occurrences of unexpected waves 
than those implied by the conventional return periods not conditioned on a reference threshold. Ultimately, 
it appears that a rogue wave that is also unexpected would have a lower occurrence frequency than that of a 
usual rogue wave. As specific applications, the Andrea and WACSIS rogue wave events are examined in detail. 
Both waves appeared without warning and their crests were nearly 2-times larger than the surrounding 0{ 10) 
wave crests, and thus unexpected. The two crest heights are nearly the same as the threshold h^ j.jpa 1.6//^. 
exceeded on average once every 0.3 • 10® waves, where Hs is the significant wave height. In confi'ast, the 
Andrea and WACSIS events, as both rogue and unexpected, would occur slightly less often and on average 
once every 3-10® and 0.6-10® waves respectively. 


1. Introduction 


A rogue wave is defined as such if the crest-to-trough 
height is at least 2.2 times the significant wave height Hs 
or if the crest height exceeds the threshold 1 .25Hs, where 
Hs = 4(7 and <7 is the standard deviation of surface eleva¬ 
tions ( [Dysthe et al.|2008^ . Evidences given for the occur¬ 
rence of such waves in nature include the Draupner and 
Andrea events. In particular, the Andrea wave was mea¬ 
sured on November 9 2007 by a LASAR system mounted 
on the Ekofisk platform in the North Sea in a water depth 
of = 74 m ( Magnusson and Donelan|2013) . The Draup¬ 
ner freak wave was measured by Statoil at a nearby plat¬ 
form in January 1995 ( |Haver|2001^ . In the last decade, the 
properties of the Draupner and Andrea waves have been 
extensively studied (Dysthe et al.|(2008); Osborne ( 1995^; 


[Magnusson and Donelan ( 2013^ ; Bitner-Gregersen et al.j 
( 2014) ; |Dias et al.| ( 2015j l and references therein). 

The Andrea wave occurred during a sea state with 
significant wave height Hs = 4a = 9.2 m, mean period 
To = 13.2 s and wavelength Lq = 220 m. The Andrea 
crest height is h = 1.63Hs = 15 m and the crest-to-trough 
height H = 2.3Hs = 21.1 m. The sea state during which 


* Corresponding author address: Georgia Institute of Technology 
Atlanta, GA 30332, USA. 

E-mail: fedele@gatech.edu 


the Draupner wave occurred had a significant wave height 
Hs = 11.9 m, mean period 7b = 13.1 s and wavelength 
Lq = 250 m. The Draupner crest height is h = 18.5 m 
(h/Hs = 1.55) and the associated crest-to-trough height 


H = 25.6 m (H/Hs = 2.15) (7Magnusson and Donelan 
(2013|l). Observations of such large extreme waves show 


that they tend to extend above the surrounding smaller 
waves without warning and thus unexpectedly. Eurther, 
both waves were twice as high as the immediately preced¬ 
ing as well as following groups of waves. In describing the 
unexpectedness of ocean waves, [Gemmrich and Garrettj 
( |2008| l define as unexpected a wave a-times larger than a 
set of one-sided (preceding) waves or two-sided (preced¬ 
ing and following) waves (see Eig. [^. Note that their def¬ 
inition of unexpectedness refers to the time interval of ap¬ 
parent calm before or during which a wave is much taller 
than the neighboring waves. Hereafter, the term ’’unex¬ 
pected” refers to a wave that is not anticipated by a casual 
observer as emphasized by jGemmrich and Garrettj ( |2010| l. 
Clearly, unexpected waves defined in this way are pre¬ 
dictable in a statistical sense as one can estimate the as¬ 
sociated return period or frequency of occurrence. 

Indeed, unexpected waves occur often with a small or 
average wave height, but they are rarely the largest waves 
in a record or r ogue waves ([Gemmrich and Garrett|2010^ . 
In this regard, jGemmrich and Garrett (20081 performed 
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Monte Carlo simulations of second order nonlinear seas 
characterized with the typical JONSWAP ocean spectrum 
and initial homogeneous random conditions. They esti¬ 
mated that a wave with height at least twice that of any of 
the preceding 30 waves occurs once every 10^ waves on 
average. Also unexpected crest heights are more proba¬ 
ble than unexpected wave heights as they occur on aver¬ 
age once every 7 • 10^ in Gaussian seas and once every 10“* 
waves in second-order nonlinear seas (see Fig. 2 in|Gemm-| 
[rich and Garrett|p008| )). Thus, their numerical predictions 
indicate that in weakly nonlinear seas unexpected waves 
occur frequently and more often than in Gaussian seas. 

Further, |Gemmrich and Garrett| ( |2008[ ) noted in their 
simulations that among the unexpected waves 2-times 
larger than the surrounding 30 waves, only about q = 
10 — 20% were rogues. With reference to second order 
crest heights, this means that in a sample population of 10^ 
waves a set of 100 waves are unexpected, as they occur 
once every 10"^ waves on average. However, only about 
10 — 20 waves of the set have crest heights that are rogue, 
i.e. larger than \3AHs as the rogue threshold adopted 
by |Gemmrich and Garrett| ( [2008[ ). This implies that unex¬ 
pected wave crests that are rogue would occur less often, 
i.e. once every 10^ waves on average. Further, the percent¬ 
age q of rogue occurrences can be interpreted as the prob¬ 
ability that the crest of an unexpected wave exceeds the 
threshold 1.34//^. Consequently, unexpected crest heights 
larger than 1 34Hs would occur rarely. 

The preceding results provide the principal motiva¬ 
tion here to consider a statistical model for describing 
unexpected waves and their rogueness. We will show 
that IGemmrich and Garrett| ( |20081 l’s definition of return 
period is unconditional. In particular, it is the harmonic 
mean of the return periods of all unexpected waves with 
any amplitude. Thus, unexpected waves of moderate am¬ 
plitude occur relatively often. However, unexpected waves 
that are rogue have a lower occurrence frequency, and this 
is in agreement with |Gemmrich and Garrett] ( |2008[ )’s nu¬ 
merical predictions. 

The remainder of the paper is structured as follows. 
First, we introduce a new theoretical model for the statis¬ 
tics of unexpected waves that accounts for both second- 
and third-order nonlinearities. We also study the effects of 
nonstationarity and stochastic dependence among succes¬ 
sive waves. In particular, we present analytical solutions 
for the return period of unexpected waves and associated 
unconditional and conditional averages for crest and wave 
heights. Then, the conceptual framework is validated by 
way of Monte Carlo simulations and the theoretical pre¬ 
dictions are compared to oceanic measurements. As a spe¬ 
cific application here, we capitalize on the numerical sim¬ 
ulations of the Andrea sea state (|Bitner-Gregersen et al.| 
|2014[ |Dias et ar]|2015^ and examine the unexpectedness 
of the Andrea wave in detail. Summary and conclusions 
follow subsequently. 


2. Statistics of unexpected waves 

Consider the exceedance probability distribution of 
wave crests characterized by third-order nonlinearities and 
described by ( [Tayfun and Fedele|2007^ 

P{x) = Pr[/! > xHs] = exp (—8vo) [l +Axq (Ax}^ — l)] , 

(1) 

where v = h/Hs is the crest amplitude h scaled by the 
significant wave height Hs = Aa and xq follows from the 
quadratic equation ( |Tayfun|1980| ) 

X = xo + 2^Xq. (2) 


Here, the wave steepness /J = A 3/3 relates to the skewness 
of surface elevations (Fedele and Tayfun 2009^ and the 
parameter 

A = A 40 -I- 2 A 22 + Ao4 (3) 

is a measure of third-order nonlinearities as a function 
of the fourth order cumulants A„,„ of the wave surface 
7] and its Hilbert transform f) (Tayfun and Fedele|2007 1 . 
|Mori and Janssen| ( |20061 l assume the following relations 
between cumulants 


A22 = A40/3, Ao 4 = A40, 


(4) 


which, to date, have been proven to hold for second-order 
narrowband waves only ( [Tayfun and Lo|[T990^ . Then, A 
in Eq. © is approximated in terms of the excess kurtosis 
A40 by 


A 


appr — 


8 A 40 


(5) 


which will be used in this work. Then, Eq. Q reduces to 


a modified Edgeworth-Rayleigh (MER) distribution (Mori 


[and Janssen|2006^. Eor realistic oceanic seas the kurtosis 


A 40 is mainly affected by bound nonlinearities (Annenkov 
[and Shrira|2014[pedele|2015b|a| |. 


Consider now a time interval ^ during which a station¬ 
ary sequence of N^, = SP jT^ consecutive waves occur on 
average. We assume that neighboring waves are stochas¬ 
tically independent. This assumption is convenient for the 
theoric development of a probabilistic model. Eurther- 
more, |Borgman| ( [r970| l argues that ”... It seems reasonable 
to assume that a wave height is at most interdependent 
with the first several wave heights occurring before and af¬ 
ter it and essentially independent with waves further back 
into the past or forward into the future”. We will show 
later that this is justified as long as the sea state is broad- 
banded so that the covariance function decays sufficiently 
rapid to zero after few wave periods and successive wave 
peaks decorrelate faster. Thus, in a sample of A(j -I-1 suc¬ 
cessive waves it is irrelevant what wave is the unexpected 
wave larger than the surrounding waves. Indeed, any wave 
in the sample could be ’’p-sided” unexpected, i.e. a-times 
larger than the previous m waves and following Na — m 
waves, with m = l,...Na/2 and p =Na/m. Eor instance. 
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Fig. 1. WACSIS measurements: the observed largest crest height is a = 2-times larger than the crests of the one-sided (two-sided) Nq ~ 50 (60) 
waves. Wave parameters — 4.16 m, T„i — 6.6 s, depth = 18 m jForristall et al.|2004^. 



Fig. 2. Unexpected crest heights in broadbanded Gaussian seas. Left panel: empirical one-sided (thin dashed line with □) and two-sided (thin 
solid line with +), Na even) unexpected wave statistics versus (solid line) predicted theoretical unconditional return period Nr in number of waves 
of a wave whose crest height is a-times larger than the surrounding Nq waves for increasing values of a = 1.5,2 and 2.5. Confidence bands are 
also shown. Right panel: empirical one-sided (thin dashed line with □) and two-sided (thin solid line with +, Na even) unexpected wave statistics 
versus theoretical predictions (solid line) of the mean crest height ha.Na ^ wave whose crest height is a-times larger than surrounding Na waves 
for a — 1.5,1.75, and 2. Sea state parameters: fully developed JONSWAP spectrum (peak enhancement factor j— 1), mean period — 8.3 
s, spectral bandwidth v = 0.35, Boccotti parameters t//* = 0.65, i/zj = 0.20 and simulated ~ 10^ waves (see left panel inset). The theoretical 
predictions accounting for the stochastic independence and dependence of successive crest heights are practically the same as the sea state is 
broadbanded. 


the last wave in the sample could be larger than the pre¬ 
ceding (one-sided) Na waves {m = Na and p = 1), or the 
central wave could extend above the preceding and follow¬ 
ing (two-sided) m=Nal2 waves (p = 2 and Na even) (see 
Fig. [TJ. Note that our definition of two-sided unexpected¬ 


ness is different than that in |Gemmrich and Garrett| ( |2^008| ) 
as they consider Na waves on each side. 

Clearly, the statistics of one- and two-sided unex¬ 
pected waves, or more generally the p-sided statistics are 
the same if stochastic independence of successive waves 
holds. On this basis, the fraction of waves n{x; ot^Na) that 
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have a dimensionless crest height h/Hs within the inter¬ 
val [x^x + dx) and that is a-times larger than any of the 
surrounding Na waves is given hy 


i{x-,a,Na)dx =\\ — P p{x)dx, ( 6 ) 


where P{x) is the exceedance prohahility given in Eq. Q 
and 


p{x) 


dP 

dx 


(7) 


is the pdf of x. Then the prohahility that the crest height ^ 
is in {x,x + dx) follows as 


Ph{x\a,Na)dx = 


nix', a,Na)dx 

n{a,Na) 


( 8 ) 


where n{a,Na) is the fraction of waves whose crest height 
is a-times larger than the surrounding Na waves, namely 

n{a,Na) = j n{x',a,Na)dx = j p{x)dx. 

° “ (9) 

By definition, the unconditional return period R or the av¬ 
erage time interval between two consecutive occurrences 
of the unexpected wave event S' is 


’ N„n{a,Na) N„n{a,Na) n{a,Na)' 

( 10 ) 

Since is the mean wave period, S occurs on average 
once every Nr waves where 


^RiO.,Na) 


1 

n{a,Na)' 


( 11 ) 


Another statistical interpretation of the unconditional 
return period Ar is as follows. Consider the average num¬ 
ber of unexpected waves nj{a,Na)Ax with a crest height 
between xj — Axf2 and Xj + Nxjl, where Ar ^ 1 is small 
and Xj are increasing amplitudes starting from xi = 0 , i.e. 
Xy+i > Xj, for j = 1, _Then, 


NR.jia,Na) = 


nj{a,Na)Ax 


is the return period of an unexpected wave whose crest 
height is nearly xj. Then, Eq. ([TT]) is approximated as 


NR{a,Na) 


1 

LJ=i«;(a,Ao)Ax 


1 

L oo _ 1 _ 

j=l N,tj{a,Na) 


, (12) 


which reveals that Nr is the harmonic mean of the return 
periods Nrj of all unexpected waves with any crest height. 

The associated mean crest height of a wave a-times 
larger than the surrounding Na waves follows from Eq. ([ 8 ]l 




= H, 


xpi,{x', a,Na)dx. 


(13) 


For comparison purposes, we also consider the standard 


statistics fimax.n, hn and hija for crest heights (Tayfun and 


Fedele 


(2007 1 ). In particular, h„ 


, is the mean maximum 


crest height of a sample of n waves 




{\-[\-P{x)Y}dx, (14) 


which admits Gumbel-type asymptotic approximations 
(Tayfun and Fedele ( 2007| i; |Fedele| ( p015a Q. Further, hn 
is the threshold exceeded by the l/« fraction of largest 
crest heights and it satisfies 


P{hnlH,) = 

n 


(15) 


where P{x) is the unconditional nonlinear probability of 
exceedance for crest heights given in Eq. 0 - The statis¬ 
tics h^jn is the conditional mean |/i > /ii/„, namely the 
average of the \/n fraction of largest crest heights 




= hn + nH, 



(16) 


One can show that hm 3 x,n is always smaller than h^/n and 
they tend to be the same as n increases (|Tayfun and Fedele| 

[20071 ). 

We also consider the standard conditional return period 
Nh ) (in number of waves) of a wave whose crest exceeds 
the threshold h = ^Hg, namely 


Nh{^) 


1 


1 


(17) 


where the exceedance probability P{^) is that in Eq.( 0 . 
From Eq. ([T3, the threshold exceeded with probabil¬ 
ity l/n implies that Nh{hn/Hs) = n, i.e. on average hn is 
exceeded once every n waves. 

Similar statistics for the crest-to-trough height y = 
PI/Hs of unexpected waves follow by replacing the crest 
exceedance probabili ty P in Eq. ([T]i with the gen eralized 
Boccotti distribution ( jAlkhalidi and Tayfun|2013| l 


PH{y)=Pr[H>yHg] = 


coexp 


4/ 


1 -I- v/-* 


1-f 


A/ 

1 -I- v/-* 




(18) 


where 

\ + \ir* 

Co = — , , 

and y/* = V^(t*) is the absolute value of the first min¬ 
imum of the normalized covariance function y/(T) = 
rj (t)ri(t + 1 ) jd^ of the zero-mean random wave process 
17 (f), which is attained at T = T* and yr* the correspond¬ 
ing second derivative ( |Boccotti|2000| i. 

The corresponding linear statistics of unexpected wave 
crests follow by setting ji = 0 and A = 0 in Eq. 0 , or 
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A = 0 in Eq. ( [TSl l for wave heights. These will hereafter he 
differentiated with the superscript L. In the following, we 
will not dwell that much on unexpected wave heights, but 
our main focus will be the statistics of unexpected crests 
in typical oceanic sea states. 

Finally, we point out that our present theory of unex¬ 
pected waves can be generalized to space-time extremes 
drawing on |Fedel^ ( |2012| i, but this is beyond the scope of 
this paper. 



Fig. 3. Conditional return period of large unexpected waves in Gaus¬ 
sian seas: (square) empirical one-sided unexpected wave statistics ver¬ 
sus (solid lines) predicted theoretical conditional return periods 
in number of waves of unexpected waves whose crest height is greater 
than ^Hs and a = 1.5-times larger than the suiTounding Na waves for 
= 0,1.0 and 1.2. Note that = 0) is the unconditional return pe¬ 
riod Nr. Sea state parameters: fully developed JONSWAP spectrum 
(peak enhancement factor y = 1), mean period T,,, = 8.3 s, spectral 
bandwidth V = 0.35, Boccotti parameters \ir* = 0.65, = 0.20 and 

simulated ~ 10® waves. The predictions accounting for the stochastic 
independence and dependence of successive crest heights are practically 
the same as the sea state is broadbanded. 


a. Stochastic dependence of successive waves 


The statistics of unexpected waves presented so far 
does not take into account the stochastic dependence of 
neighboring waves or wave groupness. Clearly, for large 
crest heights, as argued by Borgman ( | 1970) 1, one expects 
that only few neighboring crests are more or less corre¬ 
lated ( Watson) 1954[ ). To quantify this, we draw on )Fedele] 
( )2005 I and model a stationary sequence of wave crests 
jvj = as a one-step memory Markov chain. 


where each crest height xj is only stochastically dependent 
on the preceding crest height i, that is 


Since the sequence is stationary, the conditional pdf 
p{xj\xj-i) is the same for any j, say p{x 2 \x\) = 
p{xi,X 2 )/p{xi), where the crest xy precedes X 2 , and 
p{xi,X 2 ) and p{xi) are the associated joint and marginal 
pdfs. 

On these assumptions, following )Fedele) l )2005) l the frac¬ 
tion of waves n{x',a,Na) that have a dimensionless crest 
height within the interval {x,x + dx) and that is a-times 
larger than any of the surrounding Na waves is given by 


n{x;a,Na)dx = j^Pr^ 


A 

.^2 < — 

a 


X 

Xl < — 

a 


N„-l 


•Pr^V 2 < — Xl =x'j p{x)dx. 


(19) 


where 


Pr V 2 < - 
V a 


x \ _ Jo^ /o“ p{xi,X 2 )dxidx 2 
■^1 ^ ~ ) — 


a, 


and 


Pr 2^2 < — 

a 


Xl =x] = 


ff‘p(xi)dxi 


/o“ p{x,X2)dX2 


p{x) 


Then, the return period R{a,Na) of unexpected wave 
crests follows from Eqs. © and O- Clearly, if succes¬ 
sive waves were stochastically independent, p{xi,X 2 ) = 
p{xi )p{x 2 ) and Eq. reduces to Eq. © for the station¬ 
ary case. 

The theoretical probability structure of two consecutive 
wave crests is known for Gaussian processes and it is given 
by the bivariate Rayleigh distribution ( Fedele|2005)l 


Pr{xuX 2 ) = exp 


-8 


+X2 

l-k^ 


k I6k 


XiX2 


l-k^ 

( 20 ) 

where Io{y) is the modified Bessel function and the param¬ 
eter k = Xj/(T 2 ) = y /2 with the abscissa of the second 
absolute maximum of the normalized covariance func¬ 


tion V^(t) of the zero-mean random wave process (Fedele 
|2005) |. Further, the marginal pdf 

poo 

Pr{xi) = pw{xi,X2)dx2 = l(3Xitxp[-%x\) 


is the univariate Rayleigh distribution. As xj/^ tends to 
zero, successive crests become stochastically independent 
and the sea state tends to be broadbanded. Thus, we expect 
that stochastic dependence of waves is dominant in very 
narrowband sea states, where —>■ 1. In particular, our 
numerical simulations discussed later in section [^suggest 
that the dependence of consecutive crests in Gaussian seas 
is dominant when y/j >0.7. This condition corresponds 
to unrealistic oceanic sea states characterized by a Jon- 
swap spectrum with a peak enhancement factor y > 100 
and very narrowbanded as the spectral bandwidth v < 0.1. 
For typical oceanic seas, v ~ 0.3 —0.5 and ~ 0.2 —0.5, 


p{xj\xj-uXj- 2 ,---X 2 ,Xi) = p{xj\xj-i). 
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Fig. 4. Unexpected wave heights in Gaussian seas: (Left panei) predicted theoreticai unconditionai return period Nr in number of waves (soiid 
iine) versus empiricai one-sided (+) and two-sided (□, Na even) statistics as a function of the number N^, of surrounding waves for a = i.5; 
(center panei) predicted mean unexpected wave height versus observations as a function of the return period Nr. For comparison purposes, 

predicted mean wave height J/max.iVjj, conditionai mean and (right panei) threshoid //a's versus observations (circies) are aiso shown. Sea 

state parameters: fuiiy deveioped JONSWAP spectrum (peak enhancement factor 7 = 1), mean period = 8.3 s, spectrai bandwidth V = 0.35, 
Boccotti parameter Y* = 0.65 and simuiated ~ iO® waves. 


and successive waves can be assumed as stochastically in¬ 
dependent. 

Note that, the joint pdf of consecutive Gaussian wave 
crests in Eq. (|20l) can be generalized to account for second- 


order bound nonlinearities following Fedele and Tayfun 


( [2009^ , but this is beyond the scope of this work. Since 
second-order bound harmonics are phase-locked to the 
Fourier components of the linear free surface, we expect 
the classical Tayfun’s (1980) enhancement of successive 
linear crest amplitudes, but their dependence should be un¬ 
affected by second-order nonlinearities. 


b. Nonstationarity 

The statistics of unexpected waves formulated so far is 
valid for stationary sea states. In nonstationary seas, as 
those during storms, our present theory can be formalized 
as follows. From Eq. the pdf of an unexpected wave 
crest height h generalizes to 


Ph{x',a,Na\bi,.. .bm) with respect to the random vari¬ 
ables bj, that is 

Ph{x-,a,Na)Ns = Phix-,a,Na\h)'’, 

where the vector b = [b\^.. .bM] and the labeled overbar 
denotes statistical average with respect to b only. Taylor- 
expanding around the mean b = ,..., up to second 

order, yields 

Ph{x-,a,Na\h'^ + (b-b)'’+ 

Ph{x-a,Na)NS- ^ 

(b-b)^H(b) (b-b) 

( 22 ) 

where the superscript T denotes matrix transposition, the 
vector g has entries 

_ dph{x-,a,Na\b\,...bM) 

■ db, ,^5 


Ph{x;a,Na)Ns = j J Ph{x-,a,Na\bi,.. .bM) 


( 21 ) 


■p{bu...bM)dbi ...dbM: 


and the Hessian matrix 


[H(b)] „=//„ = 


d^Ph{x;a,Na\bu...bM) 


dbdu 


b=b 


where ure M time-varying wave parame¬ 

ters, e.g. a,/i,A 4 o,A 22 and A 04 , the conditional pdf 
Pk{x\a,Na\b\^.. .bM) is the stationary pdf in Eq. 
for given values of bj and p {bi,.. .bM) is the joint pdf 
of the parameters, which encodes their time variabil¬ 
ity. Eq. (|2T|) can be interpreted as the average value of 


Taking the averages in Eq. ^22) yields 

Ph{x;a,Na\bi,...bM) + 

Ph{x-,a,Na)NS - 


3^Ph(Na,Na\bi,...bi,, 


Br 


b=b 


(23) 
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where 

^rs — {pr ^r) {ps 

are correlation coefficients, in particular B^r = CT^ is the 
variance of br- These can he easily estimated from the 
nononstationary time series. Thus, p/, is the sum of i) the 
pdf in Eq. ([^ evaluated using the mean parameters b) and 
ii) additional terms that account for the spreading of the 
parameters from their mean. A similar formula can he ob¬ 
tained for n[a,Na) in Eq. The statistical moments of 
Pf, can then he obtained by integrating Eq. ( |2^ and the 
nonstationary return period follows from Eq. O- 
In our applications (see section]^, time wave measure¬ 
ments at a point are subdivided in a sequence of opti¬ 
mal 30-min intervals during which the sea state can be 
assumed as stationary. We observed that shorter time 
intervals lead to unstable estimates of higher order mo¬ 
ments, whereas longer intervals violate the stationarity 
assumption. The variability of the standard deviation a 
was taken into account by normalizing the surface height 
measurements in each 30-min interval by the respective 
observed a. In our data analysis, wave parameters are 
estimated as the average values over the available time 
record. Then, the statistics of unexpected waves can be 
based on Eq. ( |2^ , where the terms accounting for non- 
stationarity are neglected. 



Fig. 5. The role of stochastic wave dependence to the unexpected¬ 
ness of crest heights in narrowband Gaussian seas: (thin dashed line 
with □) empirical one-sided unexpected wave statistics versus (solid 
lines) predicted theoretical unconditional return periods Nr in number 
of waves for (thin line) independent and (thick line) dependent crest 
heights of a wave whose crest height is a-times larger than the sur¬ 
rounding Na waves; a — 1.5,1.75 and 2. Confidence bands are also 
shown. Sea state parameters: Gaussian spectrum with spectral band¬ 
width V = 0.1 (similar to a Jonswap spectrum with peak enhance¬ 
ment factor 7 ~ 300), mean period — 8.3 s, Boccotti parameters 
1 ^* = 0.94, y/J = 0.82 and simulated ~ 10^ waves. 


3. Are rogue waves really unexpected? 

Our interest is to describe statistically the occurrence of 
rogue waves with crest heights larger than 1 .25Hs ( |Dysthe| 



Fig. 6 . Unexpected crest heights in unidirectional second-order ran¬ 
dom seas. Empirical one-sided (thin dashed lines with □) unexpected 
wave statistics versus (thick solid lines) predicted theoretical uncondi¬ 
tional return period Nr in number of waves of a wave whose crest height 
is a-times larger than the suiTounding Na waves for increasing values 
of a = 1.5,1.75,2,2.25 and 2.5. Confidence bands are also shown. 
Sea state parameters: fully developed JONSWAP spectrum (peak en¬ 
hancement factor Y— 1), mean period T,„ — 8.3 s, spectral bandwidth 
V = 0.35, Tayfun steepness — 0.06 and simulated ~ 10^ waves. The 
theoretical predictions accounting for the stochastic independence and 
dependence of successive crest heights are practically the same as the 
sea state is broadbanded. 


|et al.|[2008| ). For example, observations indicate that the 
Andrea rogue wave appeared without warning suddently, 
attained a crest height hohs = 1.62//j, and it was as nearly 


two-times larger than the surrounding 0(30) waves (Mag- 
[nusson and Donelan|2013^ . Thus, the Andrea wave is un¬ 
expected in accordance with the definition of |Gemmrich| 
[and Garrett| ( |2008l l. However, as it will be discussed later 
in section an application of our present theory using 
Eq. {n} predicts that a wave with a crest height at least 
twice as that of any of the surrounding Na = 30 waves 
occurs on average once every Nr ~ 10^ waves. This is 
clearly observed in the left panel of Fig. Further, 
the right panel of the same Figure shows that the actual 
Andrea crest height is nearly the same as the threshold 
/iq 3.106 ~ l.6Hs exceeded by the 1 /( 0 .3- 10 ®) fraction of 
largest crests. Eq. ( fTT] ) also suggests that the Andrea wave 
is likely a rare event as the crest threshold 1.6Hs is ex¬ 
ceeded once every A/, = 0.3 • 10® waves on average. In 
contrast, our present theory predicts that the Andrea event 
would occur relatively often as an unexpected wave, i.e. 
on average once every ~ 10 "^ waves. 

The difference in occurrence rates is explained by first 
noting that the return period Nr is the average time in¬ 
terval between two consecutive waves whose crest height 
h, of any possible amplitude, is a-times larger than the 
surrounding Na wave crests. In other words, Eq. 0 re¬ 
veals that Nr is the harmonic mean of the return periods 
of all unexpected waves of any crest amplitude, and it is 
smaller than the return period of large (rare) unexpected 
waves. Thus, unexpected waves as defined by |Gemm-| 
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[rich and Garrett| ( |2008[ ) occur relatively often with small 
or moderate amplitude. However, unexpected waves that 
are rogue are rare, in agreement with their numerical pre¬ 
dictions (see also |Gemmrich and Garrett| ( |2010| i). 

To quantify the difference in occurrence frequencies 
of small and large unexpected waves, it is natural to de¬ 
fine the conditional return period NR{^,a,Na) of an un¬ 
expected wave whose crest height h exceeds the threshold 
^Hs and it is a-times larger than the surrounding Na wave 
crests. This is given by 


NR{^-a,Na) = 


1 NR{a,Na) 

J^n(x-,a,Na)dx Ph{^;a,Na)’ 


where 


Ph{x\a,Na)= I ph{s\a,Na)ds 


(24) 


(25) 


Nr{^)> 


1 


1 


p{x)dx P{^) 


= Nh{^), (26) 


lar conclusions hold for the WACSIS rogue wave (see sec¬ 
tion |^. 


4. Verification and comparisons 

a. Monte Carlo simulations of Gaussian seas 


is the exceedance probability of the unexpected crest 
height h from Eq. •HI- Clearly, for given a and Na the 
conditional return period Nr{^) is always greater than 
the unconditional Nr for any <^ > 0, and they are the 
same if <^ = 0. The left panel of Fig. 10 shows that the 


Andrea rogue wave as an unexpected wave that exceeds 
= I. 6//5 would occur rarely, i.e. on average once ev¬ 
ery Nr{^ = 1.6) ~ 6 • 10®. Instead, unexpected waves of 
any amplitude occur more often, and on average once ev¬ 
ery Nr ~ 10^. 

Clearly, the Andrea wave is both rogue and unexpected, 
i.e. its crest is larger than the crests of surrounding waves 
and it exceeds the threshold 1.25//^ ( [Dysthe et al.||2008) . 
What is the occurrence frequency of such a bivariate event 
in comparison to being only rogue as an univariate event? 

From Eq. (|24li the following inequality holds 


where we have used n(jr; a,Na) < p{x) from Eq. ©■ Here, 
Nh{^) is defined in Eq. ( [TtI i) as the standard unconditional 
return period (in number of waves) of a wave whose crest 
exceeds the threshold h = ^Hs- Thus, a wave whose crest 
is both larger than ^Hs and unexpected (as being larger 
than the surrounding waves) has a lower occurrence fre¬ 
quency than a wave whose crest is just larger than the same 
threshold. 

The preceding results imply that a rogue wave that is 
also unexpected has a lower occurrence frequency than 
just being rogue. For example, for the Andrea sea state 
the return period of a crest larger than = 1.6Hs is 
Nh{hn) = 0.3 • 10®. This is smaller than the return period 
Nr of an unexpected wave exceeding the same threshold, 
i.e. Nr{^ = 1.6) ~ 6 • 10® (see left panel of Fig. 10 1 . Simi¬ 


Drawing on Gemmrich and Garrett ( 2008| l, we per¬ 
formed Monte Carlo simulations of a Gaussian sea de¬ 
scribed by the average JONSWAP spectmm with peak en¬ 
hancement factor 7=1. The sea state is broadbanded with 
mean period r„, = 8.3 s, peak period Tp = 10 s, spectral 
bandwidth v ~ 0.35 and Boccotti parameters y/* = 0.65, 
y/2 =0.3 (see covariance function in the panel inset of 
Fig|2|- A long time series of wave surface displacements 
was randomly generated containing a total of ~ 10® waves, 
from which unexpected waves were sampled. As the sea 
state is broadbanded, our theoretical predictions can be 
based on Eqs. (|^ and ( [TT| i assuming the stochastic inde¬ 
pendence of successive crest heights. 

The left panel of Fig. shows the empirical return pe¬ 
riod Nr = RjPyn in number of waves of both one-sided 
(thin dashed line) and two-sided (thin solid line) unex¬ 
pected wave crests as a function of the surrounding Na 
waves for different values of a {Na even for the two-sided 
statistics). The two statistics are roughly the same with 
two-sided unexpected waves slightly less frequent than the 
one-sided waves. Note that for the two-sided unexpect¬ 
edness IGemmrich and Garrett| ( |20081 l consider Na waves 
on each side, thus their two-sided return period is larger 
than ours. Shown in the right panel of Fig. are also 
the empirical statistics of mean crest heights in compar¬ 
ison to our theoretical predictions for stochastically inde¬ 
pendent waves. In particular, we note that the mean crest 
height of two-sided unexpected waves is slightly smaller 
than that of one-sided waves, especially as a increases. 
Further, in Fig. there are shown the predicted condi¬ 
tional return periods Nr{^) (solid lines) of an unexpected 
wave whose crest height is greater than ^Hs for ^ = 0,1.0 
and 1.2 (a = 1.5). Note that Nr{^ = 0) is the uncondi¬ 
tional return period Nr. We find a fair agreement with the 
empirical one-sided unexpected wave statistics (squares). 
For a = 2 and Na = 30 our predicted return period is 
Nr ~ 6 ■ 10^ and i n fair agreement with the line ar predic¬ 
tions (~ 7 • lO'^) by Gemmrich and Garrett ( 2008| l as shown 
in their Fig. 2. As regard to unexpected crest-to-trough 
heights, our theoretical model fairly predicts the empirical 
wave height statistics from simulations as clearly seen in 

Fig-E 

In the above comparisons, the fair agreement with our 
theoretical predictions indicates that the stochastic inde¬ 
pendence of waves holds approximately as the sea state 
is broadbanded. However, in very narrowband seas the 
stochastic dependence of neighboring waves cannot be 
neglected. Indeed, consider a linear sea state character¬ 
ized with a Gaussian spectrum with spectral bandwidth 
V = 0.1. This is similar to an unrealistic Jonswap spectrum 
with peak enhancement factor 7 ~ 300. From the panel in¬ 
set of Fig. E the Boccotti parameters are y/* = 0.94 and 
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Fig. 7. Left panel: WACSIS, predicted theoretical nonlinear unconditional return period Nj^ in number of waves (solid line) of a wave whose 
crest height is a-times larger than the surrounding Na waves, linear predictions (dash lines) and empirical one-sided observed statistics (□) for 
a = 1.5 and 2. Confidence bands ai‘e also shown. Right panel: same for TERN measurements. Statistical parameters are taken from |TayfunH2Q06) ; 
|Tayfun and Fedel^j2QQ7^ . 


y/* = 0.81 indicating a strong correlation between con¬ 
secutive waves. Indeed, from the same figure the empir¬ 
ical one-sided (square) unexpected wave statistics tends 
to agree with our predicted theoretical return period 
for dependent waves (thick solid line) computed using 
Eqs. and Instead, our predictions for indepen¬ 
dent waves (thin solid line) are less conservative, where 
we use Eqs. (|^ and ( [TT] ). 


b. Monte Carlo simulations of second-order random seas 


Drawing on [Tayfun and Fedele (2007 1 , we performed 
Monte Carlo simulations of unidirectional second-order 
broadband random seas in deep water described by the 
same average JONSWAP spectmm introduced in the pre¬ 
vious section for simulating Gaussian seas. The associated 
Tayfun (1980) steepness /i = A 3/3 ~ 0.06, where A 3 is the 
skewness of surface elevations ( |Fedele and Tayfun|2009 1 . 
Our theoretical predictions are based on Eqs. (|9|l and (|ll]i 
and assume the stochastic independence of successive 
crest heights as the sea state is broadbanded. 

In Fig. it is shown the comparison between the em¬ 
pirical return period Nr = R/T,,, in number of waves of 
one-sided (squares) unexpected wave crests and theoreti¬ 
cal predictions from our model as a function of the sur¬ 
rounding Na waves for different values of a. For a = 2 
and Na = 30 our predicted second-order return period 

~ 2 • 10^ is shorter than the linear counterpart (~ 6 • 10^) 
for Gaussian seas (see Fig. [2)i as non linearities enhance 
crest heights (|Tayfun and Fedele|2007[|Fedele and Tayfun] 



Fig. 8 . WACSIS unexpected wave crest heights: predicted theo¬ 
retical nonlinear (solid line) and linear (dash line) mean heights ha.Na 
and ^ function of the number Na of surrounding waves ver¬ 

sus empirical one-sided statistics (squares) for a — 1.5. Horizontal line 
denotes the observed maximum crest height 1.6271.5. Wave parameters 
Hs — 4.16 m, — 6.6 s, depth d — 18 m fFomstall et aL|2004). Aver- 
age wave parameters are taken from |Tayfun| l |20Q6k|Tayfun and Fedele| 
(2007) , in particular skewness A 3 ~ 0.23 and excess kurtosis A 40 ~ 0.11. 


20091. Further, our second-order predictions fairly agree 


with those by |Gemmrich and Garrett|p008| t in their Fig. 2. 
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For example, they predict a slightly shorter nonlinear pe¬ 
riod ~ lO'* for a = 2 and Na = 30. This is because 
their second-order correction for crest heights is based on 
the narrowband assumption of the sea state. This yields a 
slightly overestimation of crest heights shortening Nr. In 
contrast, our simulated sea states are based on the exact 
second-order solution for unidirectional broadband waves 
in deep water (|Tayfun|1980^. 


c. Oceanic observations 


We will analyze two data sets. The first comprises 9 h of 
measurements gathered during a severe storm in January, 
1993 with a Marex radar from the Tern platform located in 
the northern North Sea in a water depth oid= 167 m. We 
refer to IForristaill ( |2000| i for further details on the data set, 
hereafter referred to as TERN. The second data set is from 
the Wave Crest Sensor Intercomparison Study (WACSIS) 
( |Forristall et al.| ( |2004l i). It consists of 5 h of measurements 
gathered in January, 1998 with a Baylor wave staff from 
Meetpost Noordwijk in the southern North Sea (average 
water depth d = 18 m). |Tayfun| (|2006|l and|Tayfun and| 
|Fedele| ( |2007| l elaborated both data sets and provided ac¬ 
curate estimates of statistical parameters, especially skew¬ 
ness and fourth-order cumulants which will be used in this 
work. The data analysis indicates that the statistics of un¬ 
expected waves can be based on Eq. ( |23| ), where the 
terms accounting for non-stationarity are neglected. Eur- 
ther, successive waves can be assumed as stochastically 
independent as the both sea states are broadbanded as in¬ 
dicated by their estimated covariance functions (see panel 
insets in Eig.|7]i. 

As regard to WACSIS measurements, the left panel of 
Eig.|7]compares the theoretical nonlinear return period Nr 
(solid line) of unexpected wave crests a-times larger than 
the surrounding Na waves, the respective linear predictions 
(dashed line) and the WACSIS empirical one-sided statis¬ 
tics for a = 1.5,2 (dashed line with □). The right panel 
of the same figure shows similar comparisons for TERN. 
The observed occurrence rates are close to the theoretical 
predictions, indicating that the assumption of stochastic 
independence of waves holds approximately. It is noticed 
that nonlinearities tend to reduce the return period of unex¬ 
pected waves and increase their mean crest amplitudes. In 
particular, in the left panel of Eig. |^we compare our pre¬ 
dicted nonlinear (solid line) and linear (dash line) mean 
crest heights h{a,Na) and {a,Na) versus the WACSIS 
empirical one-sided statistics (□) for a = 1.5. Clearly, our 
linear predictions underestimate the observed crest ampli¬ 
tudes, as expected. Indeed, it is well established that non- 
linearities must be accounted for to obtain reliable statis¬ 
tics of unexpected waves (Tayfun [19801 Eorristall||2000 


Tayfun and Eedele|2007[|Eedele and Tayfun 2009 |Gemm-| 


rich and Garrett 20111. Similar trend is also observed 


for the WACSIS rogue wave as evident from the center 


panel of Eig. Here, there are shown our nonlinear pre¬ 
dicted mean crest height /imax,VRj conditional mean 
and mean unexpected crest height ha=2j^a versus their lin¬ 
ear counterparts. The right panel of the same figure de¬ 
picts the nonlinear threshold in comparison to its lin¬ 
ear counterpart. The nonlinear and linear predictions for 
the Andrea rogue wave are also shown in Eig. 

We observe that the empirical statistics tend to deviate 
from the theoretical predictions for large values of a and 
Na- In particular, for both TERN and WACSIS we could 
not produce statistically stable estimates of extreme val¬ 
ues for Na > 10 when a > 1.5 due to the limited number 
of waves in the time series (0(10^) waves in comparison 
to the 10^ waves of the simulated Gaussian seas). Nev¬ 
ertheless, the agreement between our present theory and 
observations is satisfactory and it also provides evidence 
that successive waves in the samples are approximately 
stochastically independent. 

5. How rogue are unexpected waves? 

WACSIS observations indicate that the actual largest 
crest hahs is 1.62//^. Eig.[T]shows that the WACSIS rogue 
wave is also unexpected as it is a = 2-times larger than 
the surrounding Na ~ 50 waves. According to our statis¬ 
tical model such unexpected wave would occur often and 
on average once every Nr = A - 10^ waves, as seen in the 
left panel of Eig. Here, we report the theoretical predic¬ 
tions of the unconditional nonlinear return period Nr as 
a function of Na using Eqs. and ([TT|i). Eurther, from 
the center panel of Eig. it is seen that the associated 
average nonlinear unexpected crest height fi(a= 2 ,iv„= 50 ) 
is about 1.35//s and smaller than the conditional mean 
h\/Ri^ ~ 1 .5Hs, which is slightly larger than the mean max¬ 
imum crest height fimax.VR = 1 -48//^ of Ak = 4 • 10^ waves. 
Note that these average values underestimate the actual 
maximum crest amplitude hohs ~ 1 62//^ observed. In con¬ 
trast, the right panel of Eig.|^shows that hobs is nearly the 
same as the threshold Iiq 3.6 = 1 .6Hs exceeded on average 
once every A/, = 0.3 • 10® waves (see Eq. GZ))- 

We have seen that a correct statistical interpretation of 
the WACSIS rogue wave as an unexpected event requires 
considering the conditional return period Ar(^) of an un¬ 
expected wave whose crest height is larger than (see 
Eq. (|24l)). In particular, the left panel of Eig. [9] depicts 
plots of Nr{^) as a function of Na for increasing values 
of ^ = 1,1.2,1.4,1.55 and 1.6 (a = 2). Eor ^ = 1.6//s, 
we find that an unexpected wave exceeding this thresh¬ 
old and standing above Na = 50 waves would occur rarely 
and once every Nr{^ = 1.6) ~ 0.6 • 10®, in contrast to the 
smaller unconditional value Ar ~ 4 • 10^. 

In summary, the WACSIS wave crest as both unex¬ 
pected and rogue, i.e. two-times larger than Na = 50 sur¬ 
rounding waves and exceeding the 1.6Hs, would occur 
once every Nr = 0.6 • 10® waves on average. In contrast. 
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Fig. 9. WACSIS rogue wave: (Left panel) predicted nonlinear theoretical return periods in number of waves, of unexpected crest heights 

greater than and a — 2-times larger than the surrounding Nq waves for = 0,1.0,1.2,1.4,1.55 and 1.6 (solid lines) and (square) empirical 
one-sided unexpected wave statistics. Dashed vertical line denotes return period values at Na — 50. (Center panel) predicted nonlinear mean crest 
height /Zmax,iv^T conditional mean and average unexpected crest height ha^2,Nn versur their linear counterparts as a function of number of 

waves Nr. Empirical conditional mean /i]is also shown (circles). (Right panel) predicted (solid line) and empirical (circles) nonlinear threshold 
versus its linear counterpart as a function of Nr. Dashed vertical lines denote values at = 4 • 10"^,0.3 • 10^ and 0.6 • 10^. The horizontal 
line denotes the observed maximum crest height 1.62/fy. Average wave parameters are taken from |Tayfun|j2006) ; |Tayfun and Fedele|j200'7) , in 
particular skewness A 3 ~ 0.23 and excess kurtosis A 40 ~ 0.11. 



Fig. 10. Andrea rogue wave: (Left panel) predicted nonlinear theoretical return periods in number of waves, of unexpected crest heights 

greater than ^Hs and a — 2-times larger than the surrounding Na waves for = 0,1.0,1.2,1.4,1.55 and 1.6. Dashed vertical line denotes return 
period values at Na — 30. (Center panel) predicted nonlinear mean crest height /Jmax.A'/}. conditional mean and average unexpected crest 

height ha^2.Na versus their linear counterparts as a function of number of waves Nr. Empirical conditional mean /^i/a/^ is also shown (circles). 
(Right panel) predicted (solid line) and empirical (circles) nonlinear threshold versus its linear counterpart as a function of Nr. Dashed vertical 
lines denote values at Nr — 2- 10"^, 0.3 • 10^ and 3 • 10^. The horizontal line denotes the observed maximum crest height 1.63//.?. Wave parameters 
Hs — 9.2 m, Tm — 13.2 s, depth d — 70 m jMagnusson and Donelan|2013) , skewness A 3 — 0.15 and excess kurtosis A 40 — 0.1 (Dias et aL|2Q15) . 


the WACSIS wave as a rogue event has a crest height that ceeded on average once every Ni, = 0.3 • 10^ waves. Thus, 
is nearly the same as the threshold /Jq 3.196 = 1.6//s ex- the WACSIS rogue wave has a slightly greater occurrence 
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frequency than being both rogue and unexpected since 
Nh < Nr = 0.6 • 10®. This implies that the threshold 
exceeded by the 1 /Nr fraction of the largest crests is larger 
than 1 .6Hs and nearly the same as 1 .65Hs. 


6. The Andrea rogue wave and its unexpectedeness 

As a specific application of the present theoretical 
framework, the unexpected wave statistics of the 2007 An¬ 
drea rogue wave event is examined. The actual largest 
crest height hg^s is 1.63//s and nearly two-times larger 


than the surrounding 0(30) waves (see Fig. 12 in Mag- 
[nusson and Donel^ ( |2013| l). For the hindcast Andrea sea 
state, the left panel of Fig.[T0|shows the unconditional and 
conditional nonlinear return periods Nr and Nr{^) as a 
function of Ng. In particular, according to our statistical 
model, the theoretical predictions indicate that a wave with 
a crest height at least twice that of any of the surrounding 
Ng = 30 waves occurs on average once every Nr = 2- 10^ 
waves irrespective of its crest amplitude. In contrast, an 
unexpected wave whose crest height exceeds the thresh¬ 
old 1 .6Hs occurs less often since our predicted conditional 
return period = 1. 6 ) ~ 3 • 10® is greater than the un¬ 
conditional counterpart Nr = 2 ■ 10^, as seen in the left 
panel of Fig. Furthermore, the crest height 1.6Hs is 
nearly the same as the threshold Iiq 3 ,jq 6 exceeded on aver¬ 
age once every A/, = 0.3 • 1/10® waves, as indicated in the 
right panel of the same figure. Thus, the Andrea wave has 
a greater occurrence rate than being both rogue and unex¬ 
pected since Nh < Ar = 3 • 10®, and implying the larger 
threshold Hnh = l.lGHs. 


than the surrounding 0 ( 10 ) wave crests, and thus unex¬ 
pected. The two crest heights are nearly the same as the 
threshold fios.io^ ~ 1 - 6 A! exceeded on average once ev¬ 
ery 0.3 • 10® waves. In contrast, the Andrea and WACSIS 
events would occur less often as being both unexpected 
and rogue, i.e. on average once every 3 • 10® and 0.6 • 10® 
respectively. 

Finally, we point out that our statistical model for un¬ 
expected waves supports and goes beyond the analysis 
by IGemmrich and Garrett| ( |2008[ ) based on Monte Carlo 
simulations. In particular, our statistical approach can be 
used in operational wave forecast models to predict the 
unexpectedness of ocean waves. 
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